The A^-soliton solution of the Degasperis-Procesi equation 



Yoshimasa Matsuno 
Department of Applied Science, Faculty of Engineering, 
Yamaguchi University, Ube 755-8611, Japan 

Abstract 

This paper extends the resuhs of the previous paper designated I hereafter in which the 
one- and two-soliton solutions of the Degasperis-Procesi (DP) equation were obtained and 
their peakon hmit was considered. Here, we present the general iV-soliton solution of the 
DP equation and investigate its property. We show that it has a novel structure expressed 
by a parametric representation in terms of the BKP r-functions. A purely algebraic proof 
of the solution is given by establishing various identities among the r-functions. The large 
time asymptotic of the solution recovers the formula for the phase shift which was derived 
in I by a different method. Finally, the structure of the A^-soliton solution is discussed in 
comparison with that of the Camassa-Holm shallow water wave equation. 
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1. Introduction 

The Degasperis-Procesi (DP) equation 

Ut + Sk^Ux — Uxxt + 4:UUx = SUxUxx + UUxxx (1-1) 

where u = u{x,t) and k is a real constant, is a current research interest in soHton theory. 
Here, the subscripts t and x appended to u denote partial differentiation. This equation 
has been proposed as a candidate for the completely integrable nonlinear partial differ- 
ential equations by using the method of asymptotic integrability [1]. A subsequent work 
has established its complete integrability by exhibiting some features common to the in- 
tegrable system such as the existence of the Lax pair, infinity of conservation laws and 
so on [2]. Most studies so far have been devoted to the special case k = in equation 
(1.1) [2-6]. Among them, a remarkable feature is that the corresponding equation exhibits 
peakon solutions which are represented by piecewise analytic functions. On the other 
hand, the case k ^ merits a separate treatment and it was considered quite recently in 
our work, which is designated I hereafter [7]. Indeed, performing a reduction procedure 
for the soliton solutions of the Kadomtsev-Petviashvili (KP) hierarchy, we have obtained 
the explicit one- and two-soliton solutions and explored their properties. We have found 
that as opposed to the case k, — 0, solutions recover their analytic nature and behave 
like usual solitons. Nevertheless, several new features appear which are worth remarking. 
First, the soliton velocity depends nonlinear ly on its amplitude as opposed to the usual 
linear relation. Second, the interaction process of two solitons reveals a new aspect which 
was never observed in the interaction of solitons for the typical soliton equation like the 
Korteweg-de Vries (KdV) equation. In particular, the small soliton exhibits a nonnegative 
phase shift for certain range of the amplitude parameters which is in striking contrast to 
the KdV case where the small soliton always exhibits a negative phase shift irrespective 
of values of the amplitude parameters. 

The form of the two-soliton solution presented in I is very suggestive in constructing 
the general A^-soliton solution where is an arbitrary positive integer. We emphasize 
that it has a parametric representation in terms of two r-functions each of which has the 
standard Hirota form. This fact would enable us to surmise the functional form of the 
T-functions which produce the A-soliton solution. Nevertheless, the gap between the two- 
soliton and A-soliton cases is found to be so great that the construction of the A-soliton 
solution becomes a nontrivial problem which deserves further study. 

The purpose of this paper is to provide the general A"-soliton solution of the DP 
equation (1.1) in the form of the parametric representation. The solution presented here 
will be shown to exhibit several new features when compared with existing A-soliton 
solutions of typical nonlinear evolution equations like the KdV and KP equations. The 
procedure for constructing the A-soliton solution can be accomplished along the same 
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line as in the previous work I. However, we encounter several technical difficulties in the 
general AT-soliton case which arise in reducing the proof of the solution to certain algebraic 
identities. Thus, the establishment of these identities turns out to be the main task in the 
present analysis. In conclusion, we show that as in the case of the two-soliton solution, 
the general A^-soliton solution has a simple structure expressed by two fundamental r- 
functions. Remarkably, their functional form is found to be the same as the r-function 
for the A^-soliton solution of a model shallow water wave equation introduced by Hirota 
and Satsuma [8]. In section 2, we summarize the results obtained in the previous study 
for constructing soliton solutions of the DP equation. In section 3, we carry out the proof 
of the AT-soliton solution by means of a purely algebraic procedure in which the identities 
among the r-functions will play a central role. In section 4, we investigate the asymptotic 
behavior of the A'"-soliton solution for large time and recover the formula for the phase 
shift which has been derived in I by a different method. Section 5 is devoted to discussion 
where the structure of the A'"-soliton solution is compared with that of the Camassa-Holm 
(CH) shallow water wave equation. In appendix, we apply the CKP reduction to the r- 
function of the 2A^-soliton solution for the KP hierarchy and present the r-function which 
produces the A'"-soliton solution for a member of the CKP hierarchy. 

2. Summary of the previous study 

2.1. A system of equations equivalent to the DP equation 

Here, we summarize the results associated with the DP equation while focusing on an 
equivalent system of equations to the DP equation and their soliton solutions. We describe 
only the main formulas without further explanation. For their derivation and implication, 
we refer to our work I as well as a related paper [2] . 

In order to solve equation (1.1) under the boundary condition tt(±oo,t) = 0, we 
introduce the new variable r by the relation 



where the boundary condition for r is given by r(±oo, t) = k. A crucial step for simplifying 
the analysis is to define a coordinate transformation (x, t) — > (y, r) by 
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(2.1) 



and recast equation (1.1) into the conservation form 



n + {ru)a: = 



(2.2) 



dy = rdx — rudt 



(2.3a) 



dr = dt. 



(2.36) 
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Consequently, the x and t derivatives are rewritten as 



Odd 

The transformation (2.3) is well-defined provided that r > 0. Under this condition, its 
inverse transformation is obtained by solving the system of equations for x 

xt^u{y,t). (2.56) 

Here in equation (2.5) and hereafter, we use the time variable t instead of r in view of 
(2.3b). Using the transformation (2.3), equations (2.1) and (2.2) become 

u = -r(ln r)ty + - (2.6) 

n + r\y = 0. (2.7) 

Note that by virtue of (2.7) the compatibility condition Xyt — Xfy is assured for the system 
of equations (2.5) that determine the inverse mapping. 
2. 2 Lax hierarchy of the CKP equation 

To demonstrate the complete integrability of the system of equations (2.6) and (2.7), we 
substitute equation (2.6) into equation (2.7) to obtain 

vt = -\{r\ (2.8a) 
where the new dependent variable v is defined by 

Furthermore, we can eliminate the variable r from (2.8) and derive a single nonlinear 
equation for v. It reads 



'"tyyyy +{20v-^] Vtyy + SOVyVty + \ ISVyy + 16[2v - ^ ] } Vt 



1 



2/^2 



+ [^Vyyy + 32 (^2^; '^^Vy^ (^d-'vt - ^K^^ = (2.9) 
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where = dy is an integral operator. It has been pointed out that equation (2.9) 
with v+ (1/4k^) in place of v is a member of the Lax hierarchy of the Kaup-Kuperschmidt 
(KK) (or CKP) equation. Thanks to this fact, we were able to construct explicitly the 
A^-soliton solution of equation (2.9). 
2.3. Procedure for constructing N-soliton solution 

The procedure for constructing the A^-soliton solution of the DP equation consists of two 
steps, which we shall now demonstrate. First, we write the A^-soliton solution of equation 
(2.9) in the form 

V^\{y^f)yy (2.10) 

where / is obtained by means of a reduction for the 2A'"-soliton solution of the KP hier- 
archy. See appendix for detail. It may be expressed in the form of finite sum 



/= ^ exp 



N 



N 



with 



1 ^ 1 ^ 



^— ^t-2/,o),(. = l,2,...,iV) 



gA2i_l,2j-l _ 



i2{-l,2j 



1 - 

(P^ -PJ^)(9^ - Qi) 

{Pi + (Ij){(li+Pj) 

_ {Pi - Qj){Qi-Pj) 



,(z = l,2,...,iV) 

{i,j = l,2,...,N;ij^j) 
= 1,2,. ..,N) 



iPi+Pj){(li + Qj)' 

where the parameters pi and qi in (2. lid) and (2. lie) are given respectively by 



(2.11a) 

(2.116) 

(2.11c) 
(2.11d) 
(2.11e) 



ki 



1,2, 



.,N) 



(2.11/) 



ki 



1 - 



1 \ _ L2fc2 



1,2, 



.,N). 



{2.11g) 



Here, ki and yio are the wavenumber and phase of the zth soliton, respectively and the 
conditions < nki < l,{i = 1,2, ...,N) are imposed to ensure the positivity of r. Thus, 
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the AT-soliton solution of equation (2.9) is characterized completely by the 2N parameters 
ki and yio{i = 1, 2, A/"). For N = 1,2, the explicit expressions of / are written as 

/ = 1 + — e^i + e^^i ,{N=1) (2.12a) 
ai 

/ = 1 + — e^i + —6^2 + e^^i + e^^2 
ai a2 

+ ^e^i+^^ + ?^e2^^+^^ + ^e^^+2^^ + S^e^^^+^^^iN = 2) (2.126) 
aia2 02 oi 

where 

(fei-fe2)^[(fef-feife2 + fe2V-3] 

(/ci + A:2)2[(A;? + A;iA;2 + kl)n^ - 3] ^ " ^ 

_ {2kt - klkl + 2fct)^^ - 6(fcf + fcj) 

(fci + k2Y[{kl + A;iA;2 + - 3] " ^ ^ ^ 

Second, we substitute (2.10) into (2.8a) and integrate the resultant equation by y 
under the boundary condition r(±oo, t) = k and obtain the expression of r in terms of / 

r2 = -(ln f)ty + K''. (2.13) 

To represent the solution in parametric form, we must integrate the system of equations 
(2.5). It follows from equation (2.5a) subject to the boundary condition r(±oo,t) = n 
that 

x=y+r (--l.]dy + d (2.14) 



where d is an integration constant. This constant is independent of t as confirmed easily 
using equations (2.5b) and (2.7). Substituting r from (2.13) into (2.14) and performing 
the integration, we obtain the expression of the coordinate transformation x = x{y,t), 
which, combined with (2.5b), gives a parametric representation of the A'^-soliton solution 
for the DP equation. In the simplest case oi N = 1, the explicit form of the one-soliton 
solution reads 



^(a?-l) (af-i) 

Here, the integration constant d has been set to d = ln[(Q!i + l)/{ai — 1)] with ai being 
expressed in terms of ai as 

/ (2ai-l)(ai+T) " 

Oil = \ ^ 7T7 7T (2.15c) 
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and ^1 and ai are defined respectively by (2.11b) and (2.11c). 
2.4- Remark 

The most difficult technical problem throughout the present analysis is to evaluate the 
integral in (2.14) with r constructed from the A^-soliton r-function / through the relation 
(2.13). In I, this calculation has been performed only for the case of iV = 1, 2. As will be 
developed in detail in the next section, the problem mentioned here is resolved completely 
by finding a suitable expression of r in terms of two fundamental r-functions. Another 
remark is concerned with the form of the r-function /. As shown in appendix, it can 
be represented by an equivalent determinantal form which has been used in I. However, 
/ given by (2.11) is relevant in establishing the key identities (3.13) and (3.14) below 
since the r-functions gi and g2 cannot be put into the form of determinants but they 
are expressed by pfaflBans. Consequently, we can no longer rely on the various formulas 
for determinants used frequently in proving the analogous identities associated with the 
A'"-soliton solution of the CH equation [9] . 

3. A^-soliton solution 

3.1. Parametric representation 

Now, the main result in this paper can be formulated as follows: The A/"-soliton solution 
of the DP equation (1.1) can be written compactly in a parametric representation 




(3.1a) 



x{y,t) = - + In \- d. 



« 92 



(3.16) 



Here, gi and g2 are given respectively by 



N 



N 




(i<j) 



(3.2a) 



N 



N 




(i<j) 



(3.26) 



with 




(3.2c) 



— 



{Pi-Pj){Pi - (lj){Qi-Pj)iQi - Qj) 

{Pi + P3){Pi + Qj){Qi + P3){Qi + Qj) 
{ki /Cj) \{ki kikj -\- kj^K 3] 
{ki + kjY[{kj + kikj + k])K-^ - 3] ' 



,(^,J = l,2,...,iV;^^J) 



(3.2d) 
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where the phase variable of the zth sohton is defined by (2.11b). In deriving the second 
fine of (3. 2d), we have used (2.1 If) and (2.11g). Note remarkably that gi and g2 have the 
same functional form except for the phase factors ±0j. 
3.2. Identities among T-fuctions 

We verify (3.1) in a sequence of steps. To this end, we first define the following quantities 
constructed from the r- function introduced in (A.l): 



(3.3) 



/i = T2iv(|i + ln{pi/qi),i2 + ln(p2/92), C2iv + ln(p2Jv/92iv)) (3.4) 
h = r2N{ii + ln((ji/pi),|2 + In(q'2/P2), •••,6Ar + ln(ij2Ar/p2Af))- (3.5) 

An important observation here is that the following bilinear identity holds among /, /i 
and /2: 

(3.6) 



- ffty + ftfy = «Vl/2. 

We have already met with this type of identity in a previous paper dealing with the 
iV-soliton solution of the CH equation [9]. To modify (3.6) into a known identity, it is 
sufficient to replace the expressions of /, /i and /2 given above by their determinantal 
forms using the relation [10] 



J2 



27V 



27V 



Y ^^^^ + l^il^j^ij 



i=l 



= A27V det 



+ 



Pi -Qi Pi- Qj J ^ 



(3.7a) 



<i,j<2N 



where 



27V 



Ci = ii+ Y Aij,{i = 1,2,..., 2N) 
A27V is a constant factor given by 



(3.76) 



X2N = exp 



2N 



i,J = l 
(i<J) 



2N 



WiPi - Qi) 



(3.7c) 



i=l 



and 5ij is Kronecker's delta. Then, we can confirm that (3.6) reduces to the bilinear 
identity (23) of [9]. 

Let us now apply the CKP reduction defined by (A. 4) to /, /i and /2 and subsequently 
shift the phase variables as ~ = 1) 2, N) to introduce the new quantities 

/ and g 

f = f\cKp{^i^^i-lnai) (3.8) 
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g = fi\cKp{^i ^ - In Cbi) = f2\cKp{^i ^ - In tti). 
Here, g is represented by the finite sum analogous to the expression of / 

N N N 



(3.9) 



g= ^ exp 

/i,f=0,l 



J^(/^i + + ^{2{iii'i -Hi- i'i)ln tti + ^(/Xi - ^'i)ln I ^ 1 

i=l i=l i=i V Pi/ 

1 1 ^ 



(3.10) 



Note that only the difference between / and g is the third term on the right-hand side of 
(3.10). It is easy to see that both /i and /2 reduce to the expression of g given above as 
a consequence of the CKP reduction. 

It turns out from (3.6), (3.8) and (3.9) that 



- (In /) 



ty 



K 



(3.11) 



Comparing (3.11) with (2.13) and taking account of the boundary condition r(±oo,t) = k, 
we obtain the simplified expression of r in terms of / and g 



K 



9_ 



(3.12) 



Thus, r becomes a rational function of e^*(z = 1, 2, A?"). This significant expression 
would not be derived without noticing the identity (3.6) since r from (2.13) leads simply 
to an irrational function of e^*(z = 1, 2, iV). The compact expression (3.12) enables 
us to integrate (2.14) explicitly when coupled with the following set of identities among 
f,9,9i and 

f = 9192 + n{gi,yg2 - gig2,y) (3.13) 

9 = 9i92- (3.14) 

The proof of the above identities will be done by mathematical induction similar to that 
used successfully for the proof of the A'"-soliton solutions of various soliton equations within 
the framework of the bilinear formalism [11]. 
3.2.1. Proof of (3.13) 

The identity (3.13) is established by comparing the coefficient of the factor exp Yl^=i 

+ YjiLn+i 26 {0 <n < m < N) on both sides. Let Lm,n be the corresponding coefficient 
of / and Rm,n be the one for the right-hand side of (3.13). Correspondingly, the summation 
with respect to Hi and i^i must be performed under the conditions 



lj,i + i'i = 1, {i = 1, 2, n), Hi = Vi = = n + l,n + 2, m). 
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//j = = 0, (z = m + 1, m + 2, AT). 



(3.15) 



For the purpose of performing the proof effectively, it is convenient to introduce the new 
variables crj according to the relations 



1 1 
f^i = -^C^ + CTi), = ~^*)'(^= 1'2, ...,7V) 

where Uj takes either the value +1 or —1. Consequently 

jiijij^-ViVj = ^(1 + cjjCTj), = 1,2, ...,iV) 



liiUj + liiUi = -(1 - o-jCTj), = 1,2, ...,iV). 
Taking account of the relations which follow from (2.11d)-(2.11g) and (3. 2d) 

A2i-i,2i = 2 In ai, (i = 1, 2, N) 

Mi-i,2j + Mi,2j-i = lij,ihj = 1,2, ...,-/V;z7^ j) 
as well as (3.15) and (3.17), we find that 



(3.16) 



(3.17a) 
(3.176) 

(3.18a) 
(3.186) 



Lm,n = ^ exp 



n \ ^ 

- ^ In a^ + - XI (1 + ^i^i)7ii 



m m 



i=n+l J=i 



Rm,n= l + Kj^CTi/Ci 



exp 



cr = ±l 



i=l 



(3.19) 



(3.20) 



i=n+l J=i 



To simplify (3.19) further, we use the relation 

1 



exp 



{lij + CFiaj{A2i-l,23-l - A2i-l^2j)} 



1 

2 

1 h 



(1 + a-iO-j)-^ ; ; r- + - (1 - CrjCrj) 

(Pi + 9i)(9i +Pj) 2 



{Pi - - Pj) 
{Pi+Pj){qi + qj) 



(3.21a) 
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where 



A = 2VSh^l-^K''kl {i = 1, 2, .., N) 
= {2kt - kjk] + 2k^AK^ - 6{k^ + k]), = 1, 2, N) 



dij = {ki + kjf[{kf + kikj + kj)K^ - 3], {i,3 = 1, 2, AT). 



(3.216) 

(3.21c) 
(3.21d) 



The first expression of (3.21a) is a simple consequence of (2. lid), (2. lie) and (3. 2d) while 
the second expression follows by substituting (2.11f) and (2.11g) into the first one. Thus, 



becomes 



E n 

0- = ±l i,J = l 
(i<3) 



where Cm,n is a factor independent of (Tj(z = 1, 2, n) given by 

exp 



Cm. _ rj. — 



Em sr^m 
i=n+l 2^ / = \ 7ii 



(i<3) 
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Similarly, invoking the relations 



1 1 - ^(Tik^ 



(z = l,2,...,iV) 



exp 



1 



(1 + OiGj)^,^ 



1 + KCTj/cj 

{aiki - (TjkjY[{kf - aiajkikj + ^)k'^ - 3] 



{h + kjY[{kf + kikj + /ej)«2 - 3] 
(i,j = l,2,...,iV;i^j) 
which follow from (2.11c), (3.2c) and (3. 2d), we can modify Rm,n in the form 

" " 1 - -aiki 

Rm,n — Cm,n ^ ^ {1 + K ^ ] O'iki) JJ^ 



a=±l 



i=l i=l 



1 + Kaj/ci 



(3.22a) 



(3.226) 



(3.23a) 



(3.236) 



X 



Y\ {(^ih - '7jkj)'^[{ki - aiajkikj + kj)K^ - 3]. 

(i<3) 

Define P„ by the relation 



Rm,n -^m,n — C^^^i^j^, (fci , /C2, /Cn) 



where 



1 2^iki 



a=±l 



i=l 



i=l 



~\~ Kaiki 



(3.24) 



(3.25a) 
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X 



Y\ {(^ih - crjkj)'^[{ki - aiCTjkikj + /c|)«;^ - 3] 



< jj=i 

(.i<3) 



-E n 

(T = ±l i>J = l 

(i<j) 



(3.256) 



To proceed, we make an important remark. Let Pi^n be the first term on the right- 
hand side of (3.25b) multiplied by a factor nr=i(-'- ~ '^'^kf). By virtue of the summation 
with respect to ai{i = l,2,...,n), we can see that Pi^n is symmetric and even function 
of ki{i = l,2,...,n) and -Pi,nUi=±i/« = 0. Hence, Pi,n has a factor m^i{kf - I/k^) 
in view of symmetry. This implies that the first term is not a rational function but a 
polynomial. Furthermore, since the second term on the right-hand side of (3.25b) is an 
even function of /?j(z = 1, 2, n) due to the summation with respect to ai, only the even 
power of Pi contributes to the summation. This fact indicates that the second term is 
indeed a polynomial of ^^(z = 1, 2, n). As a consequence, P„ becomes a polynomial of 
ki{i = 1,2, ...,n). 

With the above remark in mind, we now prove the identity Pn = 0(n = 1,2,..., N) by 
mathematical induction. A direct calculation shows that Pi = P2 = 0. We assume that 
Pn-2 = Pn-i = 0. Then 

n 

Pn\k,=o = 2l[kf{n^k^ - 3)Pn-i{k2,k3, ...,kn) (3.26a) 

i=2 

n 

Pn\k,=±2/n = 2 HiiK^k'^ - 1)(«2/C,2 - 4)/K'']Pn-i{k2, k^, k^) (3.266) 
i=2 

n 

Pn\k,=k, = QkliK-'kl - 4) J[{k} - klf{{kt + klk} + ki)K^ - Q{kl + /C?)«2 + 9} 



i=3 



xPn-2{k3,k3,...,kn). (3.26c) 

On account of the properties (3.26) as well as symmetry and evenness, we see that Pn can 
be factored by a polynomial 



n / A \ ^ 



(i<j) 

of A;j(z = 1, 2, n) of degree 2n^ + 2n. On the other hand, it is obvious from (3.25b) that 
Pn is a polynomial of ki{i — 1, 2, n) of degree 2n(n — 1) + 1 at most, which is impossible 
except for = 0. This completes the proof. 
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3.2.2. Proof of (3.14) 

The proof of the identity (3.14) can be done in the same way. By comparing the factor 
6^p[Z]r=i + Z^I^n+i {0 < n< m < N) on both sides of (3.14), we obtain 



a=±l 

m m 



( ~~ ) ~ + 4 ^1 + '^i'^jhij - 2 Z ^^^J'^ 

i=l ^ -^^^ i=l i,j=l i,j=l 



2i-l,2j 



+ z z 

i=n+l 3=1 



Z 

<7 = ±1 



Z + 1 Z (1 + ^i^j)7y + Z Z 



i=l 



i,j = l 



i=n+l J = i 



. (3.27) 



The calculation leading to (3.25) is applied as well to modify (3.27) further. With the aid 
of the relation 

Qi 



Pi J 1 ^ K^/C? 

which is derived simply from (2.11f) and (2.11g), we can recast (3.27) into the form 
Qniki,k2, ...,kn) = ^ n [l + - 

cr=±li=l ^ ^ 



(3.28) 



X 



Y\ {(Jih - ajkjY[{kl - GiGjkikj + /c|)k^ - 3] 



(i<3) 



- E n (i + i"'*? - n 

cr=±li=l ^ ^ i,j = 'L 

(i<3) 



(3.29) 



Note that the second term on middle line in (3.29) is an even function of /?j(z = 1, 2, n) 
so that Qn becomes a polynomial of ki{i = 1, 2, n). We now prove the identity (3.29) 
by mathematical induction. The identity holds for n = 1,2, as checked easily. We assume 
that Qn-2 = Qn-i = 0. Then 



Qn\ki=o = 2Y[ki{K'^ki - 3)Qn-l{k2,k3,...,kn) 



(3.30a) 



i=2 



Qn\k,=±2/K = 6n[(«'^' - - 4)/Av2]g„_i(/c2,A;3, ...,kn) (3.306) 

1=2 

n 

Qn\k,=k, = Gklit^'^kl - 4){^kl If Hik^ - klfiikf + klki + fc^)«^ - Q{kl + kDn'' + 9} 



i=3 



13 



xQn-2{h,k3, ■■■,kn)- 



(3.30c) 



The symmetry and evenness of Qn with respect to ki{i = 1, 2, n) as weU as the prop- 
erties (3.30) imply that Qn has a factor 

(i<3) 

whose degree is 2n^ + 2n. But, as seen from (3.29), the degree of Qn is 2n^ at most and 
hence Qn must vanish identically, completing the proof. 
3.3. Proof of the N-soliton solution 

Once the identities (3.13) and (3.14) have been established, the proof of the A'"-soliton 
solution can be done straightforwardly. Indeed, if we substitute (3.12) with (3.13) and 
(3.14) into (2.14) and integrate it with respect to y, we obtain (3.1b). The expression 
(3.1a) follows immediately from (3.1b) by a simple differentiation with respect to t while 
using (2.5b) and the constancy of d. 
3.4- Remark 

The structure of the A^-soliton solution thus obtained is worth elucidating. The present 
analysis shows that the r-functions gi and g2 defined by (3.2) characterize completely the 
A'"-soliton solution of the DP equation. Another r-function / introduced in (2.11) is then 
expressed in terms of gi and g2 as indicated by the key identity (3.13). It is interesting to 
recall that gi and g2 have the same structure as the r-function for the A'"-soliton solution 
of a model equation for shallow water waves introduced by Hirota and Satsuma [8] which 
is a member of the BKP hierarchy in view of the classification theory of soliton equations 
[12, 13]. The identity (3.13) shows that the CKP r-function / can be represented by the 
BKP r-functions gi and g2. 

Although we have addressed equation (1.1) with a nonzero k, the = is an 

interesting problem to call a special attention. As already mentioned in introduction, the 
search dealing with this issue has already been done by several authors [2-6] . In particular, 
an important problem associated with nonanalytic peakon solutions [3, 5] is how to reduce 
them from analytic soliton solutions. By taking a limit — > 0, a single peakon solution 
has been reproduced from the one-soliton solution. However, the general A'"-peakon case 
is still left as an open problem even though we have a numerical evidence showing that the 
two-soliton solution with very small k approximates quite well the two-peakon solution. 
See section 4.3 of I. 

Another limit k — > oo will also deserve remarking. It has been shown that if one 
replaces the velocity of the ith. soliton by kf in (2.11b) and takes the limit k — > oo, then 
the r-function reduced from / given by (2.11a) provides the A'^-soliton solution of the 
following fifth-order KdV equation introduced by Kaup [14] which belongs to a member 
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of the CKP hierarchy [13, 15] 



, 80 o 



(3.31) 



To clarify the structure of the limiting form of the r-function /, it is suitable to use (3.2) 
and (3.13). We find in the limit of k — > oo that gi and g2 are expanded in inverse power 
of K, as 



1 / 1 

^1 ~ /i hf + [ ^ 

1 



g2 ^ h + -hf + O 



where 



h = exp 



K, 



N 



N 



(i<j) 



ii = h{y- kft + ^t' -yio],{i= 1, 2, N) 



{hi kj) (/Cj kikj + /s • ) 



{i,j = l,2,...,N;i^j). 



(3.32a) 
(3.326) 

(3.33a) 

(3.336) 
(3.33c) 



{ki -\- kj^'^^kf -\- k^kj -\r k'j^ 

Here, t' is an auxiliary time variable which may be either set to zero or absorbed into the 
phase constant j/jo in the final stage of the calculation. Substitution of (3.32) into (3.13) 
yields 

fr^f = e- 2{hht'y - ht'hy) (3.34) 

and V with / = / in (2.10) gives the A^-soliton solution of equation (3.31). Hence, a single 
T-function h characterizes completely the A^-soliton structure. See also [16-18] as for the 
construction of soliton solutions to equation (3.31) by means of the bilinear approach. 



4. Asymptotic behavoir of the A^-soliton solution 

The asymptotic form of the iV-soliton solution has been investigated in I. Here, we provide 
an alternative but more straightforward way to derive it on the basis of the asymptotic 
form of the r- functions gi and g2- To proceed, we order the magnitude of the soliton 
velocity in the {x, t) coordinate system as 



Ci > C2 > ... > Cn 



with 



3«3 



Ci = 



1 - K^kf ' 



(i = l,2,...,iV) 



(4.1a) 



(4.16) 
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and then transform to a moving reference frame with a constant vefocity Cj. We first take 
the hmit t — > — oo with the phase variable of the zth sohton being fixed. Since other 
phase variables tend to ±00 as 



^1)^2, — *■ +00, ^i+l,^i+2, ■■■,^N — *■ -00 



(4.2) 



the r-functions gi and ^2 given respectively by (3.2a) and (3.2b) are found to have the 
leading-order asymptotics 



gi ^ exp 



g2 ~ exp 



i-l 

3 = 1 

i-l 



1 + e' 



where 



i-l 



j=i 



(4.3a) 



(4.36) 



(4.3c) 



The asymptotic forms of u and x follow simply by substituting (4.3) into (3.1). They can 
be written as 



(+) 



X-Cit- Xio ~ + In Try 



i-l 
3 = 1 



(4.4a) 
(4.46) 



where Ui{^i) is the one-soliton solution 

Uii^i) = 



nkiCi sinh 



cosh + cosh 



(4.4c) 



and Xio = yio/n. We can confirm that the above expression of Ui coincides with the 
one-soliton solution given by (2.15a). Actually, by virtue of the relations 



cosh 



1 

2ai 

a,; 



sinh 0^ = ;^^(4af-l)(af-l) 





1 




1 



(4.5a) 
(4.56) 

(4.5c) 
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which come from (2.11c) and (3.2c), we can rewrite (4.4c) as 



^(a2-l)(af-i) 
cosh + 2ai- 

which is just (2.15a). 

In the Hmit of t — > +oo, on the other hand, the expressions corresponding to (4.4a) 
and (4.4b) turn out to be as 

«~«z(^z + 7l"^) (4.7a) 
x-Cit-Xio - 4^ +ln — -— y -2V(/>j+d (4.76) 

where 

N 

7i"^= E 7ii,(^=l,2,...,iV). (4.7c) 

If one observes the behavior of the A'"-sohton solution described above in the coordinate 
system at rest, the asymptotic form of the solution as t — > ±oo is represented by a 
superposition of N single solitons 

N 

« ~ + 7^^^ {t ^ -oo) (4.8a) 

i=l 
N 



i=l 



Note that if one shifts the phase as — > — X]^=i liji ^^en the argument of Ui in 

(4.8a) and (4.8b) becomes — ■* and — 7^^^^ respectively. The resulting asymptotic 
expressions of u coincide with those derived in I by a different method (see formulas (4.45) 
and (4.50) of I). It should be emphasized, however, that these asymptotic formulas are 
valid in the (y, t) coordinate system. The new feature appears when we transform back 
to the original {x, t) coordinate system, which we shall now demonstrate. Indeed, as 
t — > — oo, the center position of the ith soliton in the (y, t) coordinate is found from (4.8a) 
as = — 7^^^"* • It then follows from (4.4b) that the trajectory of the corresponding center 
position Xc in the {x, t) coordinate is described by the relation 

x,-c,t-x,o = -^+lni^^^]-2}2<Pj + d. (4.9) 
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As t — > +00, the expression corresponding to (4.9) takes the form 



X, - c,t - x,o = + ^^[yT^) (4-10) 

* ^ ^ j=i+l 



Let Aj be the phase shift of the zth sohton which is defined by the shift of the trajectory 
of the zth sohton as t — > +00 relative to its trajectory as t — > —00. This quantity is 
evaluated simply from (4.9) and (4.10), giving rise to the result 



= i^^^^ ~ ^)+'^J2^^'^Yl 0.,(^ = l,2,...,iV). (4.11) 
If we use the formula 



i-l N 



ai + 1' 



(i=l,2,...,iV) (4.12a) 



with 



/ (2a^ - l){a^ + 1) , . , . 

^ ^/(2a. + l)(a.-l)' = ^^-'^^^ 



which is derived from (2.11c) and (3.2c), formula (4.11) reproduces the corresponding 
expression presented in I (see (4.52) of I). While the first term on the right-hand side of 
(4.11) is in accordance with the formula for the phase shift arising in the context of the 
shallow water wave equation [8], the second and third terms appear as a consequence of 
the coordinate transformation (2.3). In view of the latter terms, the characteristics of the 
interaction process of solitons are found to differ substantially from those derived from 
the shallow water theory. A full explanation about this topic has already been given in I 
in the case of the two-soliton solution. 

5. Discussion 

In conclusion, it will be worthwhile to discuss the structure of the A^-soliton solution of 
the DP equation in conjunction with that of the following CH equation [19-21] 

Ut + 2k'^Ux - Uxxt + ^UUx = '2'UxUxx + UU^xx- (5.1) 

The multisoliton solutions of the CH equation with k 7^ have been constructed by 
several different methods [9, 22-27]. In particular, we have demonstrated by using an 
elementary theory of determinants that the A^-soliton solution of the CH equation exhibits 
a parametric representation similar to the A''-soliton solution (3.1) of the DP equation. 
Indeed, it can be written as [9] 

u{y,t)^ (in^-f) (5.2a) 
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y 

x{y,t) = - + ln — + d, 
« Jl 



(5.25) 



where /i and /2 arc r-functions whose structure is essentially the same as that of the r- 
functions constructing the A^-soliton solution of a model equation for shallow water waves 
introduced in the context of the inverse scattering transform (1ST) method [26]. Thus, 
the difference between the r-functions gi, (72 and /i, /2 turns out to be the most important 
issue to be explored. Although /i and /2 have the determinantal expressions [9], gi and 
g2 are expressed by a finite sum as indicated in (3.2). An important observation is that 
gi and g2 can be written in terms of pfaffians [12, 13, 29]. These r-functions satisfy the 
bilinear equation [8] for a model shallow water wave equation of Hirota and Satsuma which 
is a member of the BKP hierarchy. On the other hand, it is well known that /i and /2 
satisfy the bilinear equation [8, 25] for a model shallow water wave equation of Ablowitz 
et al which is a member of the AKP hierarchy. These inspections show clearly that the 
DP and CH equations belong to the different class of soliton equations. The discussion 
given here becomes more transparent if one compares the Lax pairs associated with both 
equations. Indeed, the spectral problem associated with the DP equation is described by 
the third order equation [2] while the one corresponding to the CH equation is the second 
order equation [20, 21]. As far as we know, however, the former spectral problem has not 
been studied sufficiently as yet. In this respect, it is an interesting problem to derive the 
A^-soliton solution of the DP equation by means of 1ST. 
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Appendix. The CKP r-function 

We introduce the KP r-function in the standard Hirota form [30] 



r2N = exp 



2N 



2N 



i=l 



■i , J — 1 

(i<j) 



(Al) 



with 



ii = ipi-qi)y + - ^) i^'^t + iio, (i = 1, 2, 2N) 
\Pi Qi J 



iA.2) 



{A.3) 



{Pi-Qj){Qi-Pj) 

where the notation X]/i=o 1 ™pliGS the summation over all possible combinations of fii = 
0,1, //2 = 0, 1, /i2iv = 0,1. While the usual definition of the r-function introduces 
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an infinite number of independent variables, only the two variables y and t are retained 
and others are set to zero in (A. 2). The CKP reduction is defined by the following 
parameterization [15] 

P2i-i = qi,P2i =Pi,hi-i = -Pi.hi = -Qi, (i = 1,2, ...,N). (AA) 

In order to obtain the A^-soliton sollution of equation (2.9), the parameters pi and qi are 
specified by (2.1 If) and (2.11g), respectively. It then follows from (A.l), (A.4), (2.1 If) 
and (2.11g) that 

l2i-i = i2i = 6, {i = 1, 2, N). {A.5) 
Here, is given by (2.11b). Note that the phase factors fjo {i = 1, 2, 2A^) have been 
rewritten as ^2i-i,o = 62,0 = {Pi + Qi)yio {i = 1,2, ...,iV). We denote the CKP r- 
function constructed above by r2N\cKP- We have already shown in I that the function 
V = (3/4)(ln f)yy with / = T2n\ckp (see (2.10)) satisfies equation (2.9). We shift the 
phase of the zth soliton as d ^ — In (or equivalently, yio — > yio + (In cii/ki)) and 
rewrite / in the form 

N 



f= exp 



1 



M=0,1 
N 



J^(//2i-l + m){^i - In Oi) 



+ - (At2i-lA*2j-1^2i-l,2j-l + l^2ilJ'2jA2i,2j] 



2 

i ,j = l 



1 ^ 

\l^2j^2i-l,2j + l^2il^2j-1^2i,2j-l) 



2 



(^.6) 



+ 2 ^(/^2i-l/^2i^2i-l,2i + /^2i/^2i-1^2i,2i-l) 

where A2i-i,2j-i{= Mi,2j) and A2i-i,2j(= ^2i,2j-i) are given respectively by (2. lid) and 
(2. lie). We replace the dummy indices /ii(z = 1,2,...,2A^) in accordance with the rule 
A*2i-i l^i,l^2i T^iii = 1,2,..., A^). If wc substitute the relations A2i-i,2i = Mi,2i-i = 
2 In tti which can be derived from (2. lie), (2.11f) and (2.11g), then we arrive at formula 
(2.11a). 

In the present analysis, we have obtained the r-function in the form of finite sum. 
We can also express it in terms of a determinant, as already demonstrated in I. To see 
the relationship between two alternative expressions of /, we use formula (3.7) with the 
phase variables shifted appropriately. In the two-soliton case, for instance, we shift 
as — In d{i = 1, 2) in (2.12b) and then multiply it by a constant factor (aia2(^)^. 

We can confirm that the resultant expression coincides with the corresponding r-function 
presented in I (see (4.18) of I). 
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